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ABSTRACT 

The origin and evolution of shear disturbances within a stratified, inviscid, incompressible 
flow are investigated numerically by a Clebsch/Weber decomposition-based scheme. In contrast 
to homogeneous flows, within which vorticity can be redistributed but not generated, the pres- 
ence of a density stratification can render an otherwise irrotational flow vortical. In this work 
a kinematic decomposition of the unsteady Euler equations separates the unsteady velocity field 
into rotational and irrotational components. The subsequent evolution of these components is 
used to investigate the influence various velocity disturbances have on both stratified and 
homogenous flows. 

In particular, the flow within a two-dimensional channel is used to investigate the evolu- 
tion of rotational disturbances, generated or convected, downstream from an unsteady inflow 
condition. Contrasting simulations of both stratified and homogeneous flows are used to dis- 
tinguish between redistributed inflow vorticity and that which is generated by a density 
stratification. 

1. INTRODUCTION 

The fluid dynamics of density stratifications, which commonly occur in atmospheric and 
oceanic flows, are of both practical and scientific interest. One significant difference between the 
behaviors of stratified and homogeneous flows is that vorticity can be generated by the presence 
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of a density gradient. More specifically, vorticity is generated when a pressure gradient exists 
perpendicular to the density stratification (Yih 1980). 

Pressure gradients are often transient and vortical flows are usually unsteady; thus the 
need to understand and investigate unsteady flows becomes apparent. Unsteady flows in gen- 
eral, and transient vortical flows in particular, are often difficult to simulate either exper- 
imentally or numerically. However, as numerical techniques become increasingly accurate and 
efficient, they are being used to investigate flow phenomena of greater complexity. This is par- 
ticularly true of vortical flows which have a significant impact on the flow around aircraft and 
through turbomachinery. 

A kinematic description is an important, yet often neglected, approach to the under- 
standing of vortical flows. The convection of material surfaces, a fundamental aspect of un- 
steady vortical flow, is rarely explored computationally. However, kinematic descriptions of 
steady vortical flows have proven fruitful, both analytically and numerically. General theories 
on the kinematics of vortical flow have been developed by Trusdell (1954) and Hawthorne 
(1966), while Yih (1960) introduced a transformation that, when coupled to Hawthorne's work, 
was used by Marris (1964) to study the generation of secondary vorticity in a stratified fluid. A 
Weber transformation was used by Goldstein (1978) to investigate weak velocity disturbances 
to the potential flow around arbitrary obstacles while Atassi and Scott (1988) extended this idea 
to the distortion of vortical waves as they convected past a series of thin airfoils. Recently the 
linking of a fluid element's displacement to its velocity change has been explored geometrically 
by Hunt and Hussain (1991), while the convection of material surfaces, an essential feature of 
any Clebsch transformation, forms the basis of Ottino's (1982) investigations into fluid mixing. 
Clearly, kinematic descriptions have proven valuable to the analytical investigation of vortical 
flow. 

Clebsch potentials have also been used to simulate flows of applied aerodynamic interest, 
where the diffusion of vorticity, coupled with numerical economy, is of paramount concern. In 
most of these calculations, a vortical distortion is added to a steady potential flow; early exam- 
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pies of this approach include Murman and Stremel's (1982) wake simulation and SteinhofF and 
Suryanarayana's (1983) convecting vortex sheet. Grossman (1983) coupled a shock-fitting 
technique to a Clebsch potential scheme and calculated supersonic flows over a number of 
conical bodies. Lacor and Hirsch (1982) and Ecer and Akay (1983) have both simulated flows 
through turbomachinery while Chang and Adamczyk (1985a, 1985b) coupled Clebsch potentials 
to the Munk-Prim substitution principle and investigated the secondary flow in a turning 
channel. Although these and many other simulations were made tractable through a Clebsch 
decomposition, their application to unsteady flows has still not been widely exploited. 

Initial conditions can dominate the evolution of an unsteady simulation and are unfortu- 
nately not a product of the Clebsch decomposition itself. This fact alone can be enough to dis- 
courage one from simulating an unsteady flow with an arbitrary Clebsch decomposition. 
However, an initial condition based on Lagrangian coordinates and the rotational component 
of the initial velocity field becomes apparent when a Clebsch decomposition is derived from the 
Weber transformation (Hunt 1987). While it may not be possible to construct a globally valid 
Clebsch decomposition for all times, the evolution of an initial transient should be simulated 
accurately. In the calculations that follow, both steady rotational and irrotational flows are 
constructed analytically, reproduced numerically, and then perturbed by a number of unsteady 
inflow velocity disturbances. 

In the present work, the incompressible Euler equations are decomposed into scale elliptic 
and hyperbolic equations while the velocity field is separated into rotational and irrotational 
components. The rotational component is constructed from a series of complex-lamellar fields 
while the irrotational component is evaluated from a scalar potential field. By decomposing the 
velocity field into rotational and irrotational components, one can identify the origins of various 
shear disturbances. By knowing the transport equations that govern these components, one can 
determine how these disturbances are convected or propagated throughout the flow. 

In each of the simulations to follow, a steady stratified flow is disrupted by velocity dis- 
turbances imposed at the inflow boundary. The interaction between these evolving velocity 
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disturbances and a statically stable density stratification is our primary interest. To allow these 
flows to evolve nonlinearly, both their rotational and irrotational velocity components are re- 
constructed during each time step of their unsteady simulation. Furthermore, the inflow density 
stratifications are held constant to avoid introducing a buoyancy drive disturbance or a 
Rayleigh-Taylor instability from the inflow boundary. Any displacement of the density inter- 
face results from an interaction with the evolving inflow velocity disturbances and is not ex- 
plicitly imposed from the inflow boundary. 

In the first series of calculations, a steady shear layer is perturbed by an oscillating inflow 
centerline. The disturbances that are generated by this forcing evolve into a series of nonlinear 
rollups that are only weakly affected by the presence of the density stratification. Consequently 
there is very little difference between the vortical structures present within the homogeneous and 
stratified flows. 

The second series of calculations simulates the vortical disturbances generated by a 
starting inflow shear. Within these simulations, a vortex is generated at the leading edge of the 
shear layer, again attributed to the redistributed inflow shear, and again similar within both the 
homogeneous and stratified fluids. While the leading edge vortex is the primary rotational 
structure within these flows, a weaker downstream shear is generated within the stratified fluid. 
The shear generated downstream is quasi-two-dimensional in appearance and is both caused 
and controlled by the initial inflow velocity disturbances. 

In the final series of calculations, an oscillating irrotational inflow is used to generate an 
unsteady shear within a stratified flow. In these simulations, an initially uniform, irrotational 
flow is forced by an oscillating irrotational inflow that renders the resulting flow both vortical 
and time-periodic. A series of spatially growing vortical structures created downstream of the 
inflow boundary simultaneously grow and diminish in time. 

In each of the simulations to follow, a Clebsch/Weber decomposition of the velocity field 
is used to identify rotational disturbances and simulate their subsequent evolution. Within a 
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stratified fluid, vorticity can be convected, redistributed, or generated throughout the flow, and 
a Clebsch decomposition is used to identify and investigate each of these effects. 


2. ANALYSIS 


The incompressible Euler equations can be written in the Cartesian coordinate system 

a s 


Dp 

Dt 


0 


( 1 ) 


duj 

dx i 


= 0 


( 2 ) 


Du i 1 dp d 

Dt P dx; dx, ^ 


(3) 


where («,v) are the Cartesian velocity components, p is density, p is pressure, and g is 
gravitational acceleration. The material derivative is written 


D_d d 
Dt dt + Uj Bxj 


(4) 


and Lagrangian coordinates ( X,Y) satisfy 


DX t 

~dT 


0 


(5) 


where the Cartesian/Lagrangian transformation matrix is defined as: 


J 


dx t 

~dXj 


( 6 ) 


Linear momentum, Eq.(3), is multiplied by the matrix J to produce the following: 
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( 7 ) 


D 1 

l dx ( ' 

\ i I 

' m \ 

1 dp 

Dt ' 

r , 

I dXj \ 

2 sy) 

P dXj 


2.1. Homogeneous Fluid 


For a homogeneous medium Eq.(l) is identically satisfied and Eq.(7) can be integrated 
as follows 



( 8 ) 


such that 



( 9 ) 


where aj is a constant of integration dependant on the material coordinates (X, F). From 
Trusdell ( i 954) and Serrin (1959) one can define the Weber transformation 


D<f> Ufa p 

~DT == ~T~ gy ~~P 


( 10 ) 


such that Eq.(9) can then be written 


dx, d<b 

L — A A 

'■ dXj J dXj 


( 11 ) 


where Aj is a constant of integration and now: 
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( 12 ) 


DAj 

~Dt 


= 0 


Eq.(l 1) can now be multiplied by the inverse of matrix / to produce the Clebsch decomposition: 


dX; <90 


Ui = Aj~- 4 


dx t dx t 


(13) 


The rotational velocity field at t = r 0 is defined as 


r0 ( 
u l 


d4> 


dxj 


t = t n 


(14) 


and if one sets X t — x t at t = t 0 , then uj = Aj, and Eqs.(12) and (13) become 


Duf 

Dt 


0 


(15) 


and 


U t =Uj 


rO 


dx t 


d4> 

dxf 


(16) 


where u' 0 is an upstream-generated shear that has been convected into the domain. Provided 
the initial condition X-, = x t is specified at t = r 0 , a homogeneous, inviscid, incompressible flow 
can now be described by the scalar convection equations, Eqs.(5) and (15); continuity, Eq.(2); 
and the Clebsch decomposition, Eq.(16). From the Clebsch decomposition, Eq.(16), vorticity 
is defined as 




( 17 ) 
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where ep is the standard permutation tensor. From Eq.(17) it becomes clear that vorticity exits 
only within the presence of the redistributed shear that was generated upstream. 

2.2. Stratified Fluid 

For a stratified medium Eq.(7) is integrated 

J (18) 

such that 


dx t 

Wj 


a ) 



J d_ 

P dXj 



(19) 


where, as in the homogeneous case, aj is a constant of integration dependant on the material 
coordinates {X, Y ). One can again define Weber transformations 


D4> um 


( 20 ) 



such that Eq.(19) can be written 


( 21 ) 


dx i _ , #1 W 

U ' 1 dXj J dXj p dXj 


( 22 ) 


where Aj is a constant of integration that satisfies Eq.(12). The inverse of matrix J multiplies 
Eq.(22) to produce the Clebsch decomposition: 
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(23) 


m ‘ = 4-^T + 


dXj dcf) i d\{/ 


dx t dx t P dx t 


A component of the rotational velocity field at t — 1 0 can be defined as 


rO_f 9<j) 1 #\ 

“ f \ U ‘ dxi P d Xi ) 




(24) 


and if one sets X t = x t at / = / 0 , then uj 0 — Aj and Eq.(23) becomes 


a dx J d4> 1 W 

U ‘ U j dx[ dx t P dxj 


(25) 


where u ^ is again an upstream-generated shear that has been converted into the domain. Pro- 
vided the initial condition X t = x t is specified at r = t 0 , a stratified, inviscid, incompressible flow 
is now be described by the scalar convection equations, Eqs.(l), (5), (15), and (20); continuity, 
Eq.(2); and the Clebsch decomposition, Eq.(25). From the Clebsch decomposition, Eq.(25), 
vorticity is defined as 




r0 

«/ 


dXj j d\ p 


dx h 


P dx h 


(26) 


which implies that vorticity can exist in the presence of the upstream-generated shear or can 
be created by the presence of the density stratification. Furthermore, as the density field is re- 
distributed throughout the flow, so too is the opportunity to generate vorticity. 


3. UNSTEADY CONVECTION SCHEME 

Convection equations (1), (5), (15), and (20) can be written as the decoupled scalar sys- 
tem: 
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H 


(27) 


dW dF dG 

dt dx dy 

This equation is solved by a cell-centered finite volume scheme that constructs a piecewise linear 
distribution over each of the finite volumes; characteristically convects information to an inter- 
mediate time level; and updates the cell-centered variables with a midpoint rule time integration. 

To insure the production of a nonoscillatory solution, at least directionally, the piecewise 
linear distribution is constructed in a nonoscillatory manner. To avoid excessive amounts of 
damping, this distribution must be at least uniformly second order accurate. A one-dimensional 
interpolation can be written 

W=Wi + S*(x- jq) (28) 

where Sf is the slope over each finite volume. The slope associated with Harten and Osher's 
(1987) uniformly second order accurate UN02 scheme is written 

Media* 0, 1/2 - W„ 

S ‘= M2 < 29 > 

where Wf + X j 2 is obtained from a nonoscillatory quadratic interpolation 


IK+ 1/2 = 0.5 (IV, + »;+,)- 0.25 D l+m 


(30) 


where 


D i+ 1/2 = minmod{W i+x - 2 W t + W t _ „ W i + 2 - 2 W l+ , + WJ (3i) 


and 


minmod(a,b) = sign(a ) max(0,s(g«(a6) min( | a | , | b | )) 


(32) 
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Surface data is evaluated at time level n + 1/2 by following characteristics back to their spatial 
locations at time level n. For positive convection speeds, the surface data at time level n + 1/2 
is written: 


ll/ n + */2 . <r x Ax 

V¥ i+\j2j "tj + °iV 2 


(' 


«/+ 


1/2.A 


Ax 


^+1/2 A 

'V 2 


(33) 


Once the surface data and their fluxes are constructed at the time level n + 1/2, the cell-centered 


values are updated by the midpoint rule 


K +l -K - A 'i 


wi+1/2 wi+1/2 

M + II2J 1 


1 / 2 ,/ 


+ 1/2 
^7 + 1/2 


+ 1/2 
«V — 1 /2 


Ax 


Ay 


(34) 


where F = ITu and G = Wv. 


4. STEADY STATE POTENTIAL SCHEME 

For a homogenous flow the continuity equation 


V • v 


d ( ro dXj d<j> 

dxj y Uy dx t dx t 


= 0 


(35) 


must be satisfied at each time step to produce the velocity potential <f>. This equation is solved 
by the approximately LU factored scheme 

[1 - M a (^x + <5y)] * [1 + + <5~)]A <t>ij = cao(8 x u + dyv)^ (36) 

where A 4>ij = <l>ij — 4>i/, H and a are scalar constants of o(l); a> is a relaxation parameter; and 
<5 + , 8~, and 8 are forward, backward, and central difference operators. For a stratified flow the 
continuity equation 
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(37) 


V • v 


d ( rO dXj d(f> 

dx; y 7 dx t dxi 


1 d\j/ 
P dxj 



must be satisfied at each time step to produce the velocity potential \j/. This equation is solved 
by the approximately LU factored scheme 


[l -Jf- (£ + . [l + -Hjj- (5- + = ao)(S x u + S y v) g (38) 

where A^» = The residuals, Eq.(35) and (37), are approximated by a finite volume 

formulation that constructs fluxes on the faces of each mesh cell. The scalar systems are then 
solved by two explicit sweeps through the domain, similar to the procedures developed by 
Anderson, Thomas, and Whitfield, (1986) and Yokota and Caughey (1988) for two factored, 
implicit time-marching schemes. The approximate LU scheme is also written within the 
framework of the multigrid method to accelerate these calculations to a steady state. 


5. ANALYTICAL STEADY SHEAR FLOW 

From Ho and Huerre (1984), a steady, quasi-two dimensional shear layer can be modelled 

u = u{\ + Rtanhp) (39) 

v = 0 (40) 

where /? = 0.5(y —y c )IQ', u = 0.5(n, + « 2 ); A u = u 1 — u x \ R = 0.5A «/«; 6 is the momentum thick- 
ness; y c is the centerline; and «, and u 2 are the velocities of the two coflowing streams. 

5.1. Homogeneous Fluid 

A Clebsch decomposition that satisfies Eqs.(2), (5), (15), (16), (39), and (40) can be writ- 
ten as 
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i/° = u (1 + R tanh ft) — u x (41) 

/° = 0 (42) 

X=x-ut(l+Rtanhp) (43) 

Y=y (44) 

<f> — u x x + — ~~ tanh /? (An + UR tanh /?) (45) 


5.2. Stratified Fluid 

For an immiscible, stratified media defined as 

P = Pi + (P 2 ~ Pi) H(y) (46) 

where H(y ) is the heaviside unit step function and p x and p 2 are the densities of the two fluids, 
a Clebsch decomposition that satisfies Eqs.(l), (2), (5), (15), (19), (25), (39), (40), and (46) can 


be written 

u' 0 = u (1 + R tanh /?) — u x (47) 

v' 0 = 0 (48) 

X = x — ut (l + R tanh P) (49) 

Y=y ( 50 ) 

r 

$ = pgyt ( 51 ) 

<j> = n,jr + n/(l + R tanh)?) • ^-y(l 4-/? tanh)?) - u,^ — gy/ (52) 
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6. NUMERICAL ASPECTS 

Shear layer calculations, characterized by « 1 =4.0m/s; u 2 = 8.0m/s; p x = 1.0/cg/m 3 ; 
p 2 = 0.9 kg/m 3 ; and 6 = 2.0mm, are performed on a 192x96 grid (Fig.l). This grid has a uniform 
streamwise spacing of 1.6mm and is algebraically stretched in the cross-stream direction. The 
cross-stream distribution has approximately 30 cells within the initial shear layer, 10 of which 
are within the momentum thickness 0. The smallest cell is 0.2mm thick and located at the 
centerline, while the computational domain is approximately 0.1m wide and 0.3m long. 

The initial conditions are constructed analytically at t = 0 and the velocity potential <f> is 
recast as a perturbation on a uniform flow. Thus 

(f) = u l x + 4> ( 53 ) 

where 0 is obtained from the continuity equation which is solved at each time step by an ap- 
proximately factored LU scheme. The LU scheme is written within the framework of the 
multigrid method to accelerate these calculations to a steady state. At each time step, a six 
multigrid level W cycle converges the average residual to 0(10" 14 ) within five iterations (Yokota 
1992). 

6.1. Homogeneous Fluid 


An inviscid noflux or parallel flow condition 


du^_ = d}_ ,0 dX 

dy dy “ dy 


0 


(54) 


is specified at both the top and bottom boundaries of the domain, while the outflow condition 


d 2 t/° _ d 2 X _ d 2 4> 
dx 2 dx 2 dy 2 



= 4> (xty b ) = 0 


(55) 
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is specified at x = x/. A steady inflow condition 


i/° = u (1 + R tanh /?) — fi x 


rc d$ n 

' - aT"° 


ax 

dx 


= 1 


( 56 ) 


can be specified at x = x 0 or modified accordingly for an unsteady simulation. 


6.2. Stratified Fluid 


The following inviscid noflux or parallel flow condition 


di/° dp d(j) rf) gx , 1 dip 

+ u — h— — — = 0 


dy dy dy 


dy P dy 


(57) 


is specified at both the top and bottom boundaries of the domain, while the outflow condition 


dj^_ _ d 2 X _ d 1 p _ d 2 4> ^ d 2 4> d f rodX_ 1 df \ 

dx 2 dx 2 dx 2 dx 2 dy 2 d y \ dy P dy J (58) 

H x tyb) == pm r 

is specified at x = jq. A steady inflow condition 


it® = u (1 + R tanh /?) — u x 


rO 


d4> d± 

dx dx 


0 


dX 

dx 


= 1 


(59) 


and Eq.(46) are specified at x = x 0 and are again modified accordingly for an unsteady calcu- 
lation. 

The analytic steady shear flows, Eqs.41-45 for a homogeneous flow and Eqs.46-52 for a 
stratified flow, have been reproduced numerically to verify the fidelity of this Clebsch/Weber 
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decomposition approach. Both flows were reproduced accurately and a detailed investigation 
of the homogeneous simulation was performed with respect to the conservation of circulation 
and enstrophy (Yokota 1991). In both of these simulations, the absence of physical viscosity 
manifests itself at the centerline of the shear layer by admitting an infinite stretching and 
shearing of the X material lines. At the centerline of the shear layer, dXjdy -» oo as t -* oo. In 
fact, it is this in viscid behavior that limits the overall accuracy of long duration or time-periodic 
calculations. On a fixed grid, the inability to resolve the infinite shearing of the X material lines 
to a given accuracy level makes the one-dimensional conditions 



+ u 


rO dX 


= 0 


for a homogeneous flow and 


(60) 


. _ 50 ' ro dX ^ 1 50 
V dy “ dy P dy 


0 


( 61 ) 


for a stratified flow increasingly difficult to satisfy at the centerline of the shear layer. These 
equations become numerically indeterminant as t-> oo. Thus the development of an invisid 
shear instability is more likely to be caused by an unbounded, numerical centerline disturbance, 
rather than by the growth of a more physically consistent initial disturbance. 


7. CENTERLINE FORCING OF A STEADY SHEAR LAYER 

To simulate the vortical rollup of a two-dimensional shear layer, without the uncertainty 
of its numerical origin, the inflow boundary conditions, Eqs.(56) and (59), are modified with an 
oscillating centerline: 


fi 


(v-£) 

29 


(62) 
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(63) 


y<=y< + K0sm(a>t) 


0) 


2nS„u 


(64) 


where S„ = 0.032 is the Strouhal number of the shear layer's natural frequency (Ho and Heurre 
1984) and AT is a scalar constant. The chosen amplitude of this centerline oscillation is one tenth 
the momentum thickness (Af=0.1) and therefore confined to within the width of the smallest 
mesh cell. The shear layer's natural frequency, 96 hz, was chosen for this oscillation, while a 
constant time step size of 0.01 ms (or 9.6 jc 10 _4 7’ where T is the period of the centerline oscil- 
lation ) was used in the time advancement. 

The initial conditions are assumed to exist for — oo < x > + oo but in Fig.2 are shown 
only within the domain encompassed by the computational grid. The initial conditions for the 
and X material fields are identical for both the homogenous and stratified fluids, while the 
density interface of the stratified flow coincides with the centerline of the shear layer. If these 
simulations were to remain steady, the u' 0 material lines and the density interface would remain 
identical to the initial conditions, while the X material field would shear continuously with the 
flow. For a steady flow, both the 4> and \j/ potential fields are initially zero and would evolve into 
the quasi-two-dimensional fields described by Eqs.51 and 52. However, in the following simu- 
lation an oscillating centerline disturbance renders the flow both two-dimensional and unsteady. 

Within both the homogeneous and stratified fuids, an initial rollup of the shear layer re- 
sults from an in viscid instability (Brown and Roshko 1974 and Saffman 1981). Given that the 
Froude number of these simulations is F r = uj{gO) 1 ^ 2 ^ 60, inertial effects will dominate those 
due to gravity and any gravitationally-induced upstream influence should be minimal. The 
vortical rollups produced within this stratified flow are not suppressed by gravity and are 
qualitatively similar to those found within its homogeneous counterpart. 

For the stratified flow, the evolution of the t/° and X material lines are shown in Figs. 3 
and 4. The evolution of the u' 0 material lines, which are equivalent to unsteady streaklines 
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passing through the inflow boundary, reveals the initial development of the small amplitude 
waves, Fig.3a-c, their subsequent steepening, Fig.3d-e, and eventual rollup, Fig.3f-h. The evo- 
lution of the density interface is comparable but not equivalent to that of the u 70 material lines 
since the inflow location of the density stratification remains constant while the centerline of the 
shear oscillates. During the development of this unsteady shear layer, the density interface un- 
dulates and steepens to produce a significant interfolding of the two fluids. 

The inviscid shearing of the X material field, Fig.4, produces a temporally growing cross- 
stream gradient that can roll up when disturbed. Within this simulation, the shear flow insta- 
bility develops from small amplitude wave to nonlinear rollup. As the X material field convects 
downstream, it both stretches and folds around the large vortical structures present within the 
flow. 

The stretching and shearing of the redistributed shear layer can be seen within the evolu- 
tion of the streamwise and cross-stream gradients of the X material field. Within the evolution 
of the streamwise gradient field, one can see the linear disturbances, Fig.5a-b, separate into a 
series of alternating regions of high and low streamwise stretching, Fig.5h. The inflow disturb- 
ances evolve into a number of spatially growing structures whose locations coincide with the 
steepening of the u ^ material lines. Since the minimum streamwise gradient is always larger 
than the specified inflow boundary condition, no streamwise contraction of the X material field 
is produced within this flow. 

The transition from small amplitude wave to nonlinear rollup is distinguishable within 
both the U 70 material field and the streamwise gradient field, while its presence is much more 
subtle within the cross-stream gradient field. The centerline of the shear layer, located at the 
maximum cross-stream gradient, bends, breaks, and ultimately evolves into a number of discrete 
structures (Fig.6c). After the initial breakage, this gradient field is drawn, pinched, and redis- 
tributed into discrete clumps (Figs.6h). These structures, formed initially at the period of the 
inflow forcing, are continually stretched and shaped as they are convected with the flow. Fur- 
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thermore, because the inflow shear dominates the evolution of this flow, the X material field 
experiences no positive streamwise shearing. 

From the evolving ^ potential field, Fig.7, one can see that disturbances are predomi- 
nantly local to the vortical structures and parallel to the density stratification. Since vorticity 
is generated only by the gradients of i// which are perpendicular to the density stratification, the 
evolution of this flow is again dominated more by the redistributed inflow shear than by the 
presence of the density stratification. Furthermore, little vorticity is generated inside the non- 
linear rollups since the i// potential field is relatively constant within their cores. 

The time histories of circulation and enstrophy over the fixed computational domain are 
global indicators of the flow's unsteadiness. When the fluid is homogeneous and the flow is 
two-dimensional, both circulation and enstrophy will remain constant over a non-material do- 
main only while the flow remains steady. The same is true for a stratified fluid only while the 
flow is both steady and one-dimensional. 

The time histories of both the stratified and homogeneous simulations are plotted with 
respect to a nondimensional time which has been normalized by the period of the inflow forcing. 
While variations do exist, both the circulation and enstrophy histories of these two flow simu- 
lations are qualitatively similar. The density stratification does not significantly influence the 
evolution of the redistributed inflow shear. The oscillating circulation histories, Fig.8, are caused 
by a cross-stream velocity component that is being generated by the inflow forcing. The slight 

differences between these circulation histories can be attributed to the oscillating location of the 

■■ ■ 

inflow shear relative to the fixed density stratification. 

The enstrophy histories, Fig.9, reflect the initial development of the small amplitude 
waves that transition into nonlinear structures after approximately two to three periods of the 
inflow forcing. Since the enstrophy histories are virtually identical, the disturbances that are 
responsible for the variations in the circulation histories have, at best, only a minor effect on 
the evolution of the redistributed inflow shear. 
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The primary vortical structures observed within this flow are similar, if not identical, to 
those seen within a homogenous flow (Yokota 1992). The qualitative behavior of these flows 
is due to the redistributed inflow shear; and while comparisons of the evolving material fields 
do not reveal any significant visual discrepancies, subtle differences do exist. Within these 
simulations, however, it is difficult to separate the evolution of the convected shear from that 
generated by the density stratification. While vorticity is being generated by the density 
stratification, its effect is weak, localized, and dominated by the redistributed inflow shear. A 
number of starting shear layers will be simulated in the following section to better investigate 
the behavior of a shear generated by a density stratification. These simulations will allow one 
to visually separate the redistributed shear from the shear generated by the density stratification. 


8. CONVECTION AND PROPAGATION OF A STARTING INFLOW SHEAR 

An inflow shear distribution is imposed upon an initially uniform flow to simulate the 
evolution of both a homogeneous and stratified starting shear layer. As this inflow shear is 
convected downstream, a vortex is generated at its leading edge. In the simulations to follow, 
the inflow shear is imposed abruptly upon the uniform flow, V</> = u x . Thus the inflow is in- 
stantaneously accelerated to the one-dimensional shear distribution described by Eqs.39 and 40. 
Since the Froude number of these flows is Fr = 60, the evolution of the leading edge vortex will 
be dominated by inertial effects, unsuppressed by gravity, and similar within both fluids. Fur- 
thermore, because the inflow density stratification is identical to that of the initial flow, no 
buoyancy-driven disturbances are explicitly introduced from the inflow boundary. 

These calculations are started with an initial X material field and density interface identi- 
cal to those specified within the previous simulations. The initial n' 0 material field is. assumed 
not to exit at t = 0 and a constant time step size of 0.01ms (or 9.6 jc 10 _4 T where T is the period 
of shear layer's natural frequency), is again used throughout these simulations. 

After 1500 time steps, the homogeneous fluid's u' 0 and X material fields are shown in 
Fig. 10, while these material fields, together with the density interface, are found in Fig. 1 1 for the 


20 



stratified fluid. The u 10 material fields, Figs. 10a and 11a, are constructed from fluid particles 
tagged with the inflow shear and reveal both its advancing front and the location of its leading 
edge vortex. The leading edge vortex evolves in a qualitatively similar manner within both the 
homogeneous and stratified fluids. 

The most significant difference between these two simulations is the presence of a down- 
stream shear found within the stratified fluid. Downstream of the leading edge vortex and along 
the density stratification, the presence of this shear can be seen within the X material field, 
Fig. 1 1 b. Furthermore, the shear's quasi-two-dimensional appearance suggests that its propa- 
gation was rapid, one-dimensional, and independant of the leading edge vortex. The starting 
shear produces an inflow disturbance that generates vorticity as it propagates along the density 
interface. This vorticity is not convected from the inflow boundary. The circulation history, 
Fig. 12, which is calculated over the fixed computational domain and normalized with respect 
to the shear layer's momentum thickness and natural frequency, gives this explanation further 
credence. While negative circulation grows linearly within both flows, an initial burst of 
vorticity is generated within the stratified fluid. Further evidence of this phenomenon can be 
obtained at a fixed downstream location by monitoring the streamwise velocity at two points 
equally spaced above and below the initial density interface. The difference between the veloc- 
ities at points located at 55% of the domain's halfwidth above and below the density interface 
are monitored at a location 90% of the domain's length downstream from the inflow boundary. 
The time history of this difference is normalized by u ave , shown in Fig. 13, and reveals an abrupt 
creation of shear within the stratified fluid. While this shear is likely to change as the leading 
edge vortex eventually approaches and passes this downstream location, it is clear that a con- 
stant shear was generated downstream of the inflow boundary. This result is in contrast to the 
homogeneous fluid within which no downstream vorticity was created. 

To gain an understanding, albeit limited, of how this downstream shear was generated, 
one can investigate the behavior of a one-dimensional velocity disturbance. Since we are prin- 
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cipally concerned with the behavior of the initial transient, a one-dimensional assumption is 
informative. 

If upon a steady one-dimensional flow, u — uq, one imposes a transient perturbation such 
that the resulting flow can be approximated as 

u = Uq + eu x + c u 2 H — (65) 

P = Po + £ Pi +£ 2 P2 + - ( 66 ) 

^ = <Ao + #1 + £ 2 <A2 + - (67) 


then after expanding this approximation within Eq.(3) and collecting all first order terms, one 
is left with the unsteady equation: 


d“i _ 1 d£i 

dt P dx 


( 68 ) 


Since the flow was initially irrotational, t/° = 0 and therefore di^o/dx = 0. By expanding Eq.(21) 
and collecting first order terms, one is left with the unsteady equation: 


#i . , 

IT +u °-eT=- p ' 


(69) 


If one defines the linearized material derivative 


D° 

Dt 



then: 


(70) 


Z)° / #i fa 
Dt \ dx ) dx 


( 71 ) 
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Thus a velocity disturbance will generate a pressure gradient that will in turn induce a change 
in the gradient of \j/ along the linearized streamwise direction. Furthermore a shear will be gen- 
erated along the density stratification since 



and the density interface is unlikely to remain undisturbed. For a stably stratified flow, where 
P 2 — P\ < 0, velocity transients can produce the following two behaviors. 

I) 



which implies that an acceleration will produce an increase in the linearized streamwise gradient 
of x and ultimately the creation of negative vorticity. This scenario corresponds to the pro- 
ceeding simulations where an inflow shear was imposed abruptly upon V(f> — u x . In those sim- 
ulations the inflow was accelerated to the shear described by Eqs.39 and 40, and negative 
vorticity was created downstream of the leading edge vortex. 


II) 


if 


du x 

~dt 


< 0 


dp\ 

dx 


> 0 


D° ( 

Dt y dx 


< 0 


which implies that a deceleration will produce a decrease in the linearized streamwise gradient 
of \j/ x and the generation of positive vorticity. This behavior is verified in the following simu- 



lation where the inflow of the initially uniform flow, V</> = u 2 , is instantaneously decelerated to 
the shear distribution described by Eqs.39 and 40. 

The initial m' 0 and X material fields, density interface, and time step size, are identical to 
those specified within the previous simulations. The u' 0 and X material fields, as well as the 
density interface, can be seen after 1500 time steps in Fig. 14. The i/° material field is con- 
structed from fluid particles tagged with the decelerated inflow shear distribution and identifies 
both the location of the advancing front and its leading edge vortex. Again the presence of the 
leading edge vortex can be seen within the X material field; however, unlike the previous simu- 
lation, the downstream-generated shear is positive in magnitude. 

For completeness, a third simulation is performed for an initially uniform flow of 
V4> = u ave . Here the upper stream is instantaneously accelerated to u 2 while the lower stream is 
simultaneously decelerated to u x . The integrated change in the inflow velocity is zero and con- 
sequently a downstream shear should not be generated. The o' 0 and X material fields, and the 
density interface, are shown after 1500 time steps in Fig. 15. The u ^ material field, which is 
constructed from fluid particles tagged with both the accelerated and decelerated inflow shear 
distribution, identifies both the location of the advancing front and its leading edge vortex. The 
development of the leading edge vortex is again similar to that found within the previous sim- 
ulations, however, there is no evidence of a downstream-generated shear. The evidence of 
vorticity. a telltale shearing of the X material field, is not found in Fig. 15b. It is also important 
to note that this behavior is independant of the stratification's placement, relative to the 
centerline of the inflow shear. Simulations in which the location of the density stratification did 
not coincide with the centerline of the inflow shear also evolved without creating a downstream 
shear. 

The circulation histories of the three stratified flow simulations are shown in Fig. 16. From 
these comparisons one can readily identify the shear that was generated when V</> = iq and 
W(j) = u 2 . This behavior is also seen within the downstream velocity differences generated within 
these simulations, Fig. 17. One can see that an equal but opposite downstream shear was gen- 
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erated when V4> = u x and V</> = w 2 . Furthermore, this shear is related to the change in velocity 
and not to a change in momentum. 

The leading edge vortex within each of these simulations displaces the density interface 
at similar downstream locations; however, the vertical depth of these displacements is most 
shallow when V<f> — u x and increases in depth from V</> = u ave to V</> = u 2 . This behavior is con- 
sistent with the increasing absolute inflow momentum change associated with these three sim- 
ulations. 

These simulations support the belief that an integrated change in the inflow velocity is 
responsible for generating the downstream shear. Buoyancy or momentum do not directly con- 
trol the creation of this shear. Moreover, when gravity was explicitly set to zero, similar results 
were produced, further suggesting that buoyancy is not the mechanism by which this down- 
stream shear is generated. 

Laminar horizontal jets in linearly stratifed fluids with Froude numbers of Fr < 1 have 
been investigated analytically by List (1971), numerically by Peyret (1976), and experimentally 
by Voropayev, Afanasyev, and Filippov (1991). Upstream effects are usually of primary interest 
within flows having these Froude numbers; these investigations, as such, focused primarily upon 
the localized influence of the jet inflow itself. Since the density stratifications within these in- 
vestigations were linear and flow visualizations were centered primarily around the inflow re- 
gion, significant evidence of a downstream-generated shear was not observed. However, despite 
the linear stratification's more subtle influence, Peyret (1976) did observe some downstream ef- 
fects: 

The fluid injected at x = 0 pushes ahead of it the ambiant fluid, which remains, to 
some extent, channelled by the effect of the buoyancy force, which tends to prevent 

downward motion Hence the perturbations created by the penetration of the jet 

are felt at very large distances ahead of it. 
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While acknowledging the dominance of buoyancy effects within low Froude number flows, it is 
still possible that this downstream phenomena was triggered by the initial inflow velocity dis- 
turbance and not by the bulk movement of fluid induced by the continuously penetrating jet. 
The preceeding simulations support the belief that a downstream shear is generated by an inflow 
velocity disturbance and not by the presence of gravity. Furthermore, the downstream shears 
were generated quicker than can be attributed to, or hindered by, buoyancy effects. 

In the absence of any redistributed inflow vorticity, the behavior of a shear generated by 
the presence of a density stratification will be investigated within the following section. A shear 
will be generated within an initially irrotational flow by a pulsed uniform inflow. Since this flow 
is free of any redistributed inflow shear, the resulting simulation should become time-periodic. 


9. SHEAR GENERATED BY AN OSCILLATING IRROTATIONAL INFLOW 

A third and final series of calculations simulates the generation of a shear by an oscillating 
irrotational inflow. Given the unsteady uniform flow 


u = (l + K sin cot) 
v — 0 


(73) 


a Clebsch decomposition that satisfies Eqs.(l), (2), (5), (15), (19), (25), (46), and (73) can be 
written: 

i/° = 0 (74) 


K> 

V 


= 0 


(75) 


' = JC — u„ ( r — COS COM 
\ ft) / 


Y=y 


(77) 
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\]/ = pxu^K sin cot + pgyt 


(78) 


4> 


u oo x-gyt + 


4K4-) 


2 rr2 

u^K 


So) 


sin 2cot 


(79) 


For the numerical simulation of this unsteady flow to remain one-dimensional, the following 
condition 


gt 


1 # 


P dy 


0 


(80) 


must be satisfied at each time step. This equation is identically satisfied if g — 0 or if p were 
constant. However, for a stratified fluid in a nonvanishing gravitational field, Eq.(80) is satisfied 
only when the density stratification, Eq.(46), is reproduced exactly. Since this stratification is 
discontinuous, it cannot be represented numerically without some diffusion and therefore dis- 
turbances are generated at this interface. The existence of gravity, coupled with an inability to 
reproduce a discontinuous density stratification, prevents this simulation from remaining one- 
dimensional. Thus, as in the case of the steady shear layer, the transition to an unsteady two- 
dimensional flow will be caused by a numerical limitation. However, one can avoid this 
numerical uncertainty by forcing the flow with an unsteady inflow condition. 

The evolution of an unsteady, irrotational inflow, Eq.(73), is simulated by assuming the 
initial conditions 


rO r0 < i a 

u = v = if/ — (p = 0 
X = x 

y=y 

<f> = U oo X + <t> 


( 81 ) 


are perturbed by the unsteady inflow condition 
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( 82 ) 


rO K) n 

u — v =— — = 0 

ox 


d<j> 

— =u 00 Ksmo)t 


where u 00 = 6 mjs, K = 0.5 is a scalar constant, p is defined by Eq.(46), and <w = 100 hz. From 
Eqs.(68)-(72) one recognizes that this unsteady inflow is capable of generating a downstream 
shear. However, unlike the previous simulations, this shear will not be quasi- two-dimensional, 
nor will it be independant of time. This simulation is the numerical counterpart to the classic 
wave-maker problem (Yih 1980, Wehausen 1991). 

For the following simulation, the initial X material field and density stratification are 
identical to those specified in the previous calculations, and a constant time step size of 
0.01 ms (or 1x10"" 3 T where T is the period of the inflow forcing) is again used throughout. 

Since the imposed inflow is irrotational, an infinite shearing of the X material lines will 
not occur, and a time-periodic behavior can be reproduced numerically. A time-periodic flow 
is achieved after approximately 8 periods of the inflow forcing and is evident within the time 
histories of circulation and enstropy. Fig. 18 and 19, respectively. Within these figures, circu- 
lation, enstropy, and time are nondimensionalized by the shear layer's natural frequency and 
momentum thickness. The circulation history is centered on a positive, nonzero mean, a be- 
havior that is not unexpected since the unsteady one-dimensional flow, Eqs.(74-79), is 
irrotational only while the density interface remains undisturbed. This asymmetric behavior 
also exists within the enstrophy history which, along with circulation, takes the approximate 
form 


r(f) = J* o)(x,y,t) da oc sin a>t + 0.2556 


(83) 
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where r(/J, E(t), and co(xj/,t) are the time dependent functions of circulation, enstrophy, and 
vorticity. These results imply that 

E{i) oc r 2 (r) (85) 

which is satisfied by the separable representation: 

(d(xy) ( sin Sr + 0.2556 ) (86) 

The behavior of this vorticity field is largely that of a spatial distribution, periodically growing 
and decaying in time. 

The evolution of this time-periodic flow can be illustrated within a sequence of results 
over one period of the inflow forcing. The evolution of the inflow (f> potential disturbances can 
be seen in Fig.20. These periodic disturbances are convected from the inflow boundary and 
distorted along the density interface. Each of these spatially growing disturbances grows and 
decays over one period of the inflow forcing. 

The unsteadiness of this flow can also be seen within the seeming erratic behavior of the 
X material field, Fig.21. Along the density interface the X material field undergoes only minor 
amounts of local shearing, both positive and negative in sign. As these spatially growing re- 
gions convect with the flow, they undergo a fairly complicated temporal pattern of growth and 
decay. First, a shear is not simply generated at the inflow boundary and convected downstream. 
Both the mechanism by which this vorticity is generated and its downstream development con- 
tinues to be controlled by the inflow velocity disturbances. Furthermore, it is clear that, relative 
to the forced shear layer, little fluid mixing occurs within this flow. No wrapping or folding of 
material lines occurs around vortical structures of any significant virulence. While vorticity is 



periodically generated in the streamwise direction, its effect on the material field is clearly 
ephemeral. 

This flow's overall structure is not immediately obvious within the X material field, which 
undergoes only, minor amounts of localized shearing and stretching. Fortunately, the same is 
not true of its gradient fields. The streamwise gradient field, Fig.22, evolves into a series of 
doublet-like structures which contain regions of both high and low streamwise stretching. The 
minimum streamwise gradient is less than the specified inflow boundary condition and implies 
that this flow, unlike the forced shear layer simulation, undergoes both streamwise contraction 
and stretching. 

The cross-stream gradient field, Fig.23, also reveals the presence of an alternating series 
of streamwise structures, positive and negative in magnitude. While these structures are con- 
vected from the inflow boundary, their individual strengths and sizes are continually fluctuating 
in response to the inflow forcing. The maximum cross-stream gradient remains constant while 
its minimum counterpart oscillates significantly. There is a clear duality to the structure of this 
flow. 

The density interface, Fig.24, oscillates in a spatially growing manner that remains linear 
in structure and is symptomatic of poor fluid mixing. This displacement, while comparable in 
amplitude to what was produced by the forced shear layer, incurs very little interfolding of the 
two fluids. While this flow is both unsteady and vortical, the mixing of the two fluids is not very 
efficient. 

T o better appreciate the relatively complicated behavior of this flow, a sequence of results 
is used to illustrate the evolution of the unsteady vortical field. The temporal locations of these 
results, relative to one period in the circulation and enstrophy histories, are shown in Figs.26 
and 27. During the time elapsed between Fig.25a-d, the flow is dominated by structures of di- 
minishing negative vorticity. Moreover, the trace structures of positive vorticity disappear, and 
only begin to reappear as negative vorticity vanishes. Fig.25e-h shows how these newly 
emerging structures of positive vorticity grow, dominate the flow, diffuse, and eventually begin 
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to disappear. It is only as they begin to vanish that their negative counterparts begin to reap- 
pear. As before, the appearance of positive vorticity is preceded by a local minimum in the 
circulation history, while local maximums are followed by the appearance of negative vorticity. 
These transitions occur near concave down-to-concave up inflection points in the enstrophy 
history where the inflow is being either accelerated or decelerated towards the mean flow. From 
these results it is clear that while positive and negative vorticity grow and decay in a manner 
diametrically opposed to one another, at no point in time is this flow irrotational. 

10. CONCLUDING REMARKS 

A kinematic decomposition of the incompressible Euler equations can be used to separate 
a velocity field into rotational and irrotational components. The rotational component is con- 
structed from a series of complex-lamellar fields while the irrotational component is evaluated 
from a scalar potential field. By decomposing the velocity field into rotational and irrotational 
components, one can identify the origins of various shear disturbances. By knowing the trans- 
port equations that govern these components, one can determine how these disturbances are 
convected or propagated throughout the flow. 

This approach was used to investigate three different flows, each of which was rendered 
both two-dimensional and unsteady by an inflow velocity disturbance. The first simulation, the 
centerline forcing of a steady shear layer, was dominated by its redistributed inflow shear and 
not significantly influenced by its density stratification. The shear layer disturbances, which 
evolve from small amplitude waves to nonlinear rollups, are qualitatively similar within both a 
homogeneous and stratified fluid. 

In the second series of calculations, a starting shear layer, significant differences between 
the homogeneous and stratified flows were found. Unique to the stratified fluid is the inflow 
velocity disturbance's ability to generate a downstream shear. Furthermore, this shear is created 
by the density stratification but is independant of gravity. 
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In the final simulation, an oscillating irrotational inflow was used to generate a stratified 
vortical flow. In the absence of any redistributed inflow vorticity, this simulation was singularly 
dependant upon the presence of the density stratification. A series of spatially growing vortical 
structures, both positive and negative in magnitude, were periodically created along the oscil- 
, lating density interface. These structures were found to periodically grow and decay at the fre- 
quency of the inflow forcing and in a manner diametrically opposed to one another. 

The strength of this kinematic approach lies in its ability to visualize fluid mixing 
quantitatively. In each of the preceding simulations, this aspect of the Clebsch/ Weber decom- 
position was useful in investigating the flows more fully. 
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Figure 1 . — Grid. 



HIN= -0. 9992 MAX= -0.6 909- INC - 0.0193 
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(a) X material field. 
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(b) U* material field. 
Figure 2. — Initial conditions. 


35 







M1N= 0.5000 HAX= 3.5000 INC= 

URO CONTOURS 

(d) t = 2.40. 

Figure 3. — U 10 material field: forced stratified shear layer. 
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(h) t = 4.32. 

Figure 3.— Concluded. 
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(d) t = 2.40. 

Figure 4.— X material field: forced stratified shear layer. 
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(h) t = 4.32. 

Figure 4. — Concluded. 
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(a) t = 0.96. 
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(d) t = 2.40. 

Figure 5. — Streamwise gradient of the X material field: forced stratified shear layer. 
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Figure 5. — Concluded. 
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(a) t = 0.96. 
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(d) t = 2.40. 

Figure 6. — Cross-stream gradient of the X material field: forced stratified shear layer. 
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Figure 6.— Concluded. 
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MIN = -0.11 IE-01 MAX - 0. 7S2E-02 INC - 0.177E-02 
PS I CONTOURS 

(c) t = 1 .92. 



MIN = - 0.10fE-01 MAX = 0.120E-01 INC = 0.213E-02 

PSI CONTOURS 


(d) t = 2.40. 

Figure 7. — ^ potential field: forced stratified shear layer. 
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MIN = -0.16+E-01 MAX - 0.113E-D1 INC = 0.263E-02 
PS I CONTOURS 

(e) t = 2.88. 



MIN * -O.J99E-OI MAX - 0.165E-01 INC - 0.298E-02 
PSI CONTOURS 

•(f) t = 3.36. 



MIN = -Om 21 9E-QI MAX - O.lfBE-OI INC - 0.39-8E-02 
PSI CONTOURS 

(g) t = 3.84. 



MIN = -0.2Q1E-01 MAX * 0.202E-01 INC = 0.383E-02' 
PSI CONTOURS 

(h) t = 4.32. 

Figure 7.— Concluded. 
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Figure 8,— Circulation history: forced shear layers. 




(a) U* material field. 



(b) X material field. 

Figure 10— Homogeneous starting shear layer - u at t = 1 .44. 
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NIN= 0.5000 MAX- 3.5000 INC = 0.5000 

URO CONTOURS 


(a) material field. 



MIN » - 1.1195 MAX= - 0.7786 INC » 0.0212 

X CONTOURS 


(b) X material field. 



(c) Density Interface. 

Figure 1 1 . — Stratified starting shear layer -V ( j ) =u 1 att = 1 .44. 
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MIN= -3.5000 MAX = -0.5000 INC= 0.5000 

URO CONTOURS 


(a) material field. 



MIN* -1.1190 MAX* - 0.7792 INC* 0.0212 

X CONTOURS 


(b) X material field. 



(c) Density interface. 

Figure 1 4. — Stratified starting shear layer -V^’s^attsl .44. 
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M1N= -1.7500 MAX= 1.7500 1NC= 0.5000 

URO CONTOURS 


(a) U 10 material field. 



MIN= -1.1107 
X CONTOURS 

MAX = -0.7005 

INC* 

0.0211 


(b)X material field. 





(c) Density interface. 

Figure 1 5.— Stratified starting shear layer - = u ave at t = 1 .44. 
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CIRCULATION 
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TIME TIME (MILLISECONDS) 

Figure 1 6.— Circulation history: stratified starting shear layers. . Figure 1 7.— Streamwise velocity difference: stratified starting 

shear layers. 




Figure 1 8.— Circulation history: oscillating irrotationai inflow. 


Figure 1 9. — Enstrophy history: oscillating irrotationai inflow. 







MIN* -1.5+18 MAX * - 1 . 2+08 INC* 0 . 0188 

PHI CONTOURS 


(a) 5 1 = 45°. 



MIN= - / . 5758 MAX * -1 . 2560 INC* 0 . 

/>/// CONTOURS 


(b) tot = 135°. 



MIN* -1.599+ MAX * - 1.2768 INC * 0.0202 

PHI CONTOURS 


(C) tot = 225°. 



MIN* -1.616+ MAX * - 1.3018 INC = 0.0197 

PHI CONTOURS 


(d) tot = 315°. 

Figure 20. — <t> potential field: oscillating irrotationa! inflow. 


52 




M I N= -1.6+51 MAX * -1.3306 INC* 0.019 ? 

X CONTOURS 

(a) W t = 45°. 



MIN- - 1.6596 
X CONTOURS 

MAX - -1.3+89 

INC* 

0.D1 9+ 


(b) cot = 135°. 





MIN = -1.6/01 

X CONTOURS 

MAX* -1.3590 

INC * 

0.019+ 


(c) <ot = 225°. 





MI N= -1.6X56 MAX- -1.3619 INC* 0.0196 

X CONTOURS 


(d) <ot = 315°. 

Figure 21 — <t> materia! field: oscillating irrotational inflow. 
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MIN = 0.9S5E 00 MAX = 0.800E 01 INC - 0.191E 01 

STREAMNISE GRADIENT OF X { 

Figure 22. — Streamwise gradient of the X material field at cot = 135°: oscillating irrotational inflow. 



MIN * -0.902E-01 MAX = 0.682E-01 INC - 0.317E-01 
CROSS-STREAM GRADIENT OF X 

Figure 23.-Cross-stream gradient of the X material field at 5 1 = 1 35°: oscillating irrotational inflow. 



Figure 24. — Density interface at ©t = 135°: oscillating irrotational inflow. 
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MIN = -0.+25E 02 MAX = 0.000 INC = 0.H-25E 01 

VORTICITY CONTOURS 

(b) wt = 63°. 



MIN - -0.22VE 02 MAX - 0.000 INC = 0.22H-E 01 

VORTICITY CONTOURS 

(C) cot = 99°. / 



MIN = -Om 32 IE 01 MAX = 0.893E 01 INC = 0.126E 01 
VORTICITY CONTOURS 


(d) cot = 126°. 

Figure 25.— Vorticity: oscillating irrotational inflow. 
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MIN = 0-000 MAX * 0.323E 02 INC - 0.323E 01 

VORTIClTY CONTOURS 

(e) wt = 171°. 



MIN = 0.000 MAX - 0.713E 02 INC = 0.713E 01 

VORTIClTY CONTOURS 

(f) cot = 216°. 



MIN = 0.000 MAX = O.SfSE 02 INC = O.SH-SE 02 

VORTIClTY CONTOURS 

(g) cot = 279°. 



MIN = -0.265E 02 MAX = 0.S06E 01 INC - 0.326E 01 
VORTIClTY CONTOURS 

(h) wt = 333°. 

Figure 25. — Concluded. 
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DEGREE 


Figure 26. — One period of the circulation history. 



DEGREE 

Figure 27. — One period of the enstrophy history. 
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